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Abstract 


Seismic shaking is an attractive mechanism to explain the destabilisation of 


regolith sl opes and the regolith mi grati on 


(Richardson et ah. 

2004; 

Mivarnoto et ah. 

2007) 


'onnd on the surfaces of asteroids 


20071 ). Here, we use a continuum 


mechanics method to simulate the seismic wave propagation in an asteroid. 
Assuming that asteroids can be described by a cohesive core surrounded by 
a thin non-cohesive regolith layer, our numerical simulations of vibrations 
induced by micro-meteoroids suggest that the surface peak ground accelera¬ 
tions induced by micro-meteoroid impacts may have been previously under¬ 
estimated. Our lower bound estimate of vertical accelerations induced by 
seismic waves is about 50 times larger than previous estimates. It suggests 
that impact events triggering seismic activity are more frequent than previ¬ 
ously assumed for asteroids in the kilometric and sub-kilometric size range. 
The regolith lofting is also estimated by a first order ballistic approximation. 
Vertical displacements are small, but lofting times are long compared to the 
duration of the seismic signals. The regolith movement has a non-linear de¬ 
pendence on the distance to the impact source which is induced by the type 
of seismic wave generating the first movement. The implications of regolith 
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concentration in lows of surface acceleration potential are also discussed. We 
suggest that the resulting surface thermal inertia variations of small fast ro¬ 
tators may induce an increased sensitivity of these objects to the Yarkovsky 

effect. 
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1. Introduction 


In-situ observations of asteroids have highlig 


ited the complexi t ies and 


varia tions in the small body surface environments flCheng et al. 


19971 : 


Kawaguchi et al 


20031). In addition to finding each of these bodies to be regolith-covered, there 
is strong evidence that this regolith is very complex and active. Evidence 


of rego 


ith motion has been observed on several as 


slides (|Thomas et al. 


Sierks et al. 


20 


Fuiiwara et al 


2001 


2002 


Veverka et al. 


lal) . crater erasure flVeverka et al. 


2001a. 


Cheng et al 


2006 1. dust “ponds” ( Veverka et al. 


eroids, includin g land- 


Rob 


2nmb 


nson et al. 


Salto et al. 


2001b 


20021 ) and particle size segregation (lYano et all. 2006J). 


2002a; 


2006 


Robinson et ah 


As impact-induced seismic shaking may trigger g. 
cesses in the dry, vacuum, microgravity environment 


obal-scale granu 


Miyamoto et al. 


ar pro 


20071), 


it is, therefore, an attractive mechanism to explain the destabilisation of 
regolith slopes and the regolith migration found on the surfaces of aster¬ 
oids. The idea of global seismic activity resulting from large impacts and 
the consequ ences o n regolith motion is not a new concept. For example, 
As phau g and Meloshl (119931 1 use hydrocodes to study the seismic resurfacing 
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ofregolith at the site of large hyper-velocity impacts. Greenber g et alj (11994 


19961 1 make reference to the idea of global jolting resulting from large impacts 


and suggest that this explains asteroid morphologies a nd cr ater populations. 


In addition to studying high-energy impacts, 


Housenl (119811 1 consider the as¬ 


teroid surface evolution produced by lower energy, non-catastrophic impact 
events. However, the focus is on the local consequences of the cratering pro¬ 
cess (e.g., burial by ejecta or excavation) rather that the global consequences 
of a small impact. 

Here, rather than looking at the global consequences of high-energy im¬ 
pacts or the behaviour of regolith at the site of the impacts, we consider the 
motion of the non-cohesive regolith layer across the entire surface as a result 


of small micro-meteroid impacts. Several models have al ready estima 


seismic activity resultin g from meteoroid impacts (e.g. ICheng et al. 


;ed the 


2002 


Richardson et al. 


2004) • In this study, we first present a new method for 


simulating the seismic wave propagation in an asteroid and then we consider 
the consequences of micro-meteoroid impacts on the regolith mobility. 

In the following, we will assume that small asteroids are composed of a 
cohesive, or monolithic, core surrounded by a layer of non cohesive, or loose, 
regolith material. We do not speculate on the type of cohesive forces at work 
in the cohesive core, or if this core is in fact composed of piles of blocks of var¬ 
ious sizes, but simply assume that the seismic waves considered in this study 
do not induce irreversible deformations of the core. In addition, only small 
amplitude seismic waves created by micro-meteoroid impacts are studied, in 
order to avoid irreversible deformations of the core. With such assumptions, 
the propagation of seismic waves in the core can be described by continuum 
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mechanics. This assumption does not prevent diffraction of seismic waves by 
heterogeneities, voids or fractures. 

The surrounding non cohesive layer is defined as a layer of loose material 
experiencing plastic (irreversible) deformations under the action of the seis¬ 
mic waves considered here. In particular, this layer is assumed to be able 
to detach from the cohesive core and subsequently experience lofting when 
the vertical accelerations induced by seismic waves exceed the local verti¬ 
cal surface acceleration. The migration of loose regolith materials has al¬ 


ready. been observed on aster oi ds sue 


Miyamoto et al. 


2007t 


i as Itokawa, Eros and Lut etia (e.g. 


Robinson et al. 


2002b 


Sierks et al. 


201lbh 


This study presents continuum mechanics seismic wave propagation simula¬ 
tions, and suggests that the surface vertical accelerations induced by micro¬ 
meteoroid impacts may have been under-estimated. As a consequence, we 
argue that impact-induced regolith motion may be more common than pre¬ 
viously thought. We quantify the regolith lofting by presenting first order 
estimates of vertical displacements and lofting time based on a ballistic ap¬ 
proximation. Finally, we consider the implications of regolith concentrations 
in surface acceleration potential lows. 


2. Seismic modelling 

2.1. Objective and previous studies 

The objective of our modelling is to provide an estimate of peak ground 
accelerations and velocities induced by direct waves for small impacts and 
various asteroid models in the sub-kilometric size range. The challenge in 
this estimate is the proper modelling of high frequency signals, because for 


4 


























small impactors (mass < 100 g) and small asteroids (diameter < 1 km) both 
the excitation at the source and the elastic response of the medium present 
dominant frequencies larger than 20 Hz. Examples of power spectral densi¬ 
ties of acceleration impulse responses for different asteroid sizes and different 
seismic wave attenuations are presented in Figure [Q (a). For asteroid sizes 
smaller than 1 km diameter, even at low quality factors (high attenuation) 
the impulse response is dominated by frequencies larger than 20 Hz. This is 
due to the fact that propagation distances are small. 

Most of the computations of seismic shaking effects on asteroids have been 


performed up to now through empi rical or ana. 


large amplitude terrestrial explo sions flCheng et al 
t hroug h hydrocode simulations (Richar ds on et al. 


ytica l relations scaled to 


2002 


2004; 


Asphaua. 


2008), or 


Walker and Huebner, 


2004al) . The majority of the studies were focused on the irreversible ground 


deformatio ns due to large impactors (lAsphaug and Meloshl . 


1994 


1993 ; 


Greenberg et al. 


19961 ). 


Two teams performed 


rydrocodes simulations of impacts on asteroid diame¬ 


ters as small as 1 km (Richardson et al 


Walker et al 


20 061) 


Richardson et al. 


2005 


Wal ke r and Huebnerl . 


2004al: 


(120051 ) suggests that small impactors 


may generate peak ground accelerations larger than local gravity on small 
asteroids. However, their hydrocode simulations are limited by their high 
frequency resolution limit at least above 50 Hz, and comparison of syn- 


thetic displacement recorc 


of 


Richardson et ah 


s with theroretical Green functions (see Fig. 8 


(120051 )) suggest that the computed asteroid responses 


are lo w passed versi o ns o: 


20 Hz. 


Walker et al 


theoretical one with a corner frequency around 


(120061 ) performed hydrocode simulations up to 800 Hz, 
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wit h a proper quantification of an explosive source flWalker and Huebnerl . 
20 04 a). However, the attenuation model and the amount of high frequencies 
removed by the hydrocodes limitations are not fully described. Moreover, 
the study was focused on the Eros case which falls outside our size range. 


Previous simulations of seismic wave propagation in a cohesive 


continuum mechanics seismic wave equations flMartin et al. 


2008 


objec 


Blitz 


with 


2009) 


predict, for objects of diameter inferior to 20 km, vertical accelerations 
at least two order of magnitudes larger than the local gravity for micro¬ 
meteoroid impacts (10 grams at 6 km/s). These simulations were performed 
by a spectral element numerical method including wave scattering by a net¬ 
work of fractures. Outputs of these models show that, despite the strong 
scattering of the wave field, the peak vertical ground acceleration is due to 
the coherent direct waves. However, the model complexity and large com¬ 
putation times make it difficult to explore various asteroid sizes and shapes 
with such simulations. That is why we present here a proper modelling of 
the high frequency elastic source signal and asteroid response through seismic 
normal modes computation. 

2.2. modelling tools and hypothesis 

Here, simulations of surface accelerations induced by micro- meteoro id im¬ 


pacts are performed by normal mode computations flZhao and Dahlenl . 


1993 


19951 1 on spherical bodies of various radii. These computations do not repro¬ 
duce the scattered wave field, but do allow an estimation of the maximum 
peak ground acceleration due to the direct waves. Because the propagation 
distances are small, and assuming that the material remains cohesive, the 
scattered wave field will generate peak ground accelerations smaller than the 
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(Larose et ah. 

2005; 

Martin et al.. 

00 

o 

o 

CM 

Blitz, 

o 

o 

CM 


estimates of peak ground accelerations will be valid if the energy loss of direct 
waves due to scattering is properly taken into account by our modelling of 
seismic attenuation. 

A simple homogeneous internal structure model is tested with characteristics 


correspondin g to Mo on mega- 


400 m/s ((Cooper et ah 


1974 


regolith: p = 


Horvath et al 


50 0 kp/m 3 , Vp = 900 m/s , Vs = 


1980 


Gagnepain-Beyneix et al. 


2006|). This model can be viewed as a worst case in terms of accelerations 


because changing from spherical to elliptical shape would focus the seismic 
energy in some particular places, and adding a low velocity layer at the 
surface would trap the seismic waves in this layer, increasing the surface ac¬ 
celerations. 

The source is described by an impact of a lOg projectile at 6 km/s at vertical 
incidence. Previous experimental and numerical studies of such small hyper¬ 
velocity impacts favour an equivalent seismic moment tensor close to the one 


of an explosion (IWalker and Huebner . 


2004a). The total seismic moment for 


sue 


by 


i an impact is estimated to be 1.8e+6 N.m from numerical experiments 


Walker and Huebnerl (j2004al) in the energy range considered here. Thi s 


value is about 10 times smaller than the est imate from 
and 10 times larger than the estimate from lCollins et al. 


Edwards et al 


(12008 ). 


( 20051) . These large 


variations are due to two poorly known parameters: the seismic amplification 
factor and the seismic efficiency. The seismic amplification factor (ratio of 
the incident impulse to the impulse transfered to the target) can v ary from 


0 .5 t o 20 depending mainly on the properties of the target (IMcGarr et al. 


19691 ). The seismic efficiency (part of the source energy converted into seismic 
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waves) varies from 0 6 to 10 2 depending o 


erties (jCpllins et ah 


2005 


Kedar et al. 


the impactor and ta rget p rop 


20121) . and impact energy (IShishkinl . 


2007 ). In our simulatio ns, t he seismic moment is set to the lowest estimated 


value of 3.77e+5 N.m (Colli ns et al. 


2005). 


Another important parameter is the equivalent seismic source duration of 
the impact, because it quantifies the time range in which the energy is con¬ 
centrated. The shorter the source duration, the lar ger the peak gro un d ac - 
celerations. Following scaling relations presented by Lognonne et al. ( 20091) . 
and using a conservative seismic efficiency of 10“ 5 , we predict a short source 
duration of 8 ms, concentrating the energy in a short time period. 

Finally, the frequency content of the signal is strongly influenced by the seis¬ 
mic attenuation. As the interiors of asteroids are cold, the instrinsic wave 
attenuation is expected to be very low, and the attenuation of direct waves 
mainly due to losses by scattering, by analogy with lunar observations. In 
addition, recent impact experiments suggest a low q uality f actor at high fre¬ 


quencies in the loose regolith layer at the source (IKedar et ah 


20121). In 


order to be conservative in our estimates we tested low quality factors of 100 
and 20 to mimic the scattering attenuation of the direct waves. 


2.3. modelling results 

Figure [0(b) presents the envelopes of surface accelerations at 90° epi- 
central distance to a micro-meteoroid impact (half way point between the 
impact and the antipode). At low quality factors, the surface accelerations 
remain during a long time and pass above the maximum surface gravity 
during tens of seconds. Decreasing the quality factor to values as low as 
20 reduces strongly the duration of the shaking. However, the energy lost 





























in scattering will arrive later, as ob served for lunar impacts ([Latham et al. . 


1970; 


Duennebier and Sutton, 


1974)). In any case, even with a high attenua¬ 


tion, the first direct wave still generates peak ground accelerations above the 
gravitational acceleration. 

In order to explore various sizes of asteroids, the computations were re¬ 
peated for asteroid diameters from 100m to 1km. The peak ground acceler¬ 
ations and velocities are presented in figure [2l Our computations up to 90° 
epicentral distance for all asteroid models overlap. It demonstrates that up 
to this distance the antipodal refocusing of the waves can be neglected, and 
it cross validates our normal modes computations for the various asteroid 
models. Except for the largest asteroids and lowest quality factor, the peak 
ground accelerations exceed the maximum surface gravity over the largest 
part of the asteroid surface. However, the peak ground velocities do not ex¬ 
ceed the escape velocity, demonstrating that these seismic vibrations generate 
regolith movements, but do not allow ejection of material. As observed on 
figure H peak ground accelerations and velocities vary more with distance to 
the source and with quality factor than with asteroid size. However, even for 
such simple internal structure models, it is not possible to provide a simple 
scaling law as a function of distance and quality factor. 


2-4- Comparison with previous estimates 

Figure [21(a) presents comparisons of estimated peak gr oun d acc elera¬ 


tions between ou r resu lts and previous ones by ICheng et al. 


Richardson et al 


(120051 ). The simple formula presented by 


(2002) and 


Cheng et al. 


(120021 ) was scaled on large amplitude terrestrial explosions, and is valid only 


in the near field range. However, it validates the fact that our peak accelera- 
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ti ons a re lower bound estimates. The estimates referring to 


Richardson et al. 


020051 ) are based on an average of seismic energy over the whole asteroid, 
that’s why they do not depend on distance, but depend on asteroid size. 
Their computations assume no attenuation, a seismic effic iency o f 10~ 5 and 


a do minant frequency of 30 Hz (equations (2) and (5) of 


Richardson et al. 


(j2005|)). Similar computations were performed with the seismic diffusion 
model (equation (14) of the same study), but due to small asteroid sizes com¬ 
pared to the minimum diffusion length estimate (125 m), the results are al¬ 
most identical. Due to the averaging over the asteroid volume the peak accel¬ 
eration values are lower than our estimates, in particular at short distances. 
Consequently we find that the minimum impactor size ne cessary to trigger 


Richardson et al. 


regolit h movement is even smaller than the estimate of 
020051 1. 

The micro-meteoroid events investigated here have a mean period be¬ 
tween imp acts sm aller than 1 yea r/km 2 for impactors more massive than 


10 grams ([Richardson et al. 


20051 ). Due to a factor ^50 increase of peak 


groun d acce l erat ion estimates over most of the asteroid surface compared 
0200511 study, impact events triggering seismic activity 


to 


Richardson et al 


are, therefore, at least 10 2 times more frequent than assumed during previ¬ 
ous studies. Despite the large error bars on the impact rate of objects more 
massive than 10 grams, kilometric scale asteroids are expected to experience 
several of these impacts during one orbit. 


2.5. Source scaling 

The peak ground accelerations and velocities obtained in this study scale 
linearly with seismic moment. The seismic moment is almost proportional 
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to the impactor kinetic energy (iWalker and Huebnerl . 2 004bl ). However, the 
peak accelerations also depend on the source duration because short sources 
of small impactors concentrate the energy in a short duration, and gener¬ 
ate high accele rations at high frequencies. Scaling relations presented by 


Log nonne et al. 


( 20091 ) suggest that the source duration scales as m3 with m 


being the impactor mass. So, the results presented here can be extrapolated 
with confidence to within a factor of 10 of the asteroid mass. 


In the case of explosive sources, the seismic moment is 


source durations are much shorter than for impacts Walker and Huebnerl . 


airly well known, bu 


2004bl h So, the extrapolation of the curves presented in Figure [2] to explo¬ 
sive sources is not possible. 

2.6. Limitations and uncertainties 

Our estimates of peak ground accelerations were constructed in order to 
be a lower bound. However, in addition to uncertainties related to the quan¬ 
tification of the seismic source, strong variations are expected due to the 
internal structure model of the asteroids. The amount of seismic scattering 
is a key parameter, because if the seismic scattering is too high the direct 
waves quantified in this study may not generate the highest peak acceler¬ 
ations. The shape of the asteroid will a’ 


and major fractures may block the waves ( Walker et al 


so stronglv 

ocus/defocus waves 

(Walker et al. 

2006; 

Martin et al. 


20081 ). These two effects create strong variations of peak ground accelera¬ 
tions on the surface. Eventually, regions with low cohesion may experience 
irreversible deformations induced by seismic waves, breaking the continuum 
mechanics hypothesis implicit to our computations. This effect is not signif¬ 
icant at low impactor energy, but is predominant at high impactor energy 
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Direct waves vertical acceleration at 90A° distance 


Apophis Model (270m diameter) 




-Q=1000 

-Q=100 

-Q=20 

-Maximum surface gravity 


✓ , v 0 5 10 15 

j K ) Time (in s) 

Figure 1: On the left, power spectral densities of acceleration impulse responses of direct 
body and surface waves at 90° epicentral distance to an impact of a lOg object at 6 km/s 
for homogeneous asteroid models of varying diameters of 100 m (red thick lines), 270 m 
(blue lines) and 1 km (black thin line). Computations are performed through normal 
mode summations as described below for different values of the quality factor of 1000 
(plain lines), 100 (dashed lines) and 20 (dotted lines). On the right, envelope of surface 
accelerations for an object of 270 m diameter for different quality factors (1000, 100, 20). 


( AsphaugL 2008!). 


3. Regolith mobility 

The resulting motion of the non-cohesive regolith layer due to the seismic 
activity may involve different types of behaviour. The first of these is a 
ballistic behaviour when the regolith particles are temporarily lofted from 
the surface. Lofting occurs when the vertical acceleration of the ground is 
larger than the ambient gravity (the gravitational attraction of the asteroid 


on a gra in and t 


Scheeres et al. 


le inertial effects that arise due to the rotation of the asteroid; 


20101 ). However, the regolith layer will only be truly lofted 


if its initial vertical velocity at the instant of ejection (given by the vertical 
velocity of the ground at the instant of ejection) is sufficient for the particle 
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Figure 2: Simulations of peak ground accelerations (a) and peak ground velocities (b) as 
a function of distance (in m) to an impact of a lOg object at 6 km/s for homogeneous 
asteroid models of following diameters: 100 m (magenta lines), 200 m (blue lines), 400 
m (green lines) and 1 km (black lines); and following quality factors: 1000 (plain lines), 
100 (dashed lines) and 20 (dotted lines). In panel (a), the maximum surface gravity 
of the various asteroid models (thick lines labelled by a steroid diameters) , and previous 


estimates presented by 1 Che ng et ajJ 120021) (CH02) and 


Richardson et al 


(1200511 (RS05) 


are also presented. In panel (b), the escape velocity of the various asteroid models (thick 
lines labelled by asteroid diameters) are also presented, assuming that the asteroids are 
non-rotating. 
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Figure 3: Vertical (on the left) and horizontal (on the right) seismic signals computed by 
our simulations. The example shown here is the seismic signal measured at a distance of 
17.5 m away from the micro-meteoroid impact on the surface of an asteroid of diameter 
100 m, where Q is assumed to be equal to 20. The records start at the arrival of the p-wave 
(the body wave), followed by the arrival of the first Rayleigh (surface) wave approximately 
20 ms later. 

to have a larger vertical displacement than the ground. Once lofted, we 
assume that the regolith layer follows a ballistic trajectory until returning 
into contact with the ground unless, of course, the initial vertical velocity is 
larger than the escape velocity in which case the regolith will be lost (note 
that in these simulations this never occurs; see Fig. [2] [b]). The regolith 
layer may be lofted several times for a given seismic event and, due to the 
low-gravity environment, even when the particles are lofted to very small 
heights, the period of free fall can be long. 

Shown in figure [3] is the resulting ground displacement, velocity and ac¬ 
celeration in both the vertical and the horizontal (radial) directions for an 
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example seismic simulation. The arrival of the p-wave (body wave) and the 
Rayleigh (surface) wave can be seen in the figure. For the range of seismic 
simulations presented above we find three different cases for the ballistic re- 
golith motion: (1) the regolith layer is lofted by the p-wave and does not 
return into contact with the ground until after the Rayleigh wave has passed 
(Fig. [4] [a]), (2) the regolith layer is lofted by the p-wave and returns into 
contact with the ground before or as the Rayleigh wave is passing and is 
subsequently re-lofted (Fig. 0] [b]), and (3) the regolith layer is never lofted. 

As the Rayleigh wave has a much higher amplitude than the p-wave, 
the particles lofted by this wave will reach larger heights and will spend 
much more time in free fall than those lofted by the p-wave. This leads to a 
counterintuitive result that, close to the micro-meteoroid impact, the regolith 
layer may be lofted less than further from the micro-meteoroid impact. This 
can be seen in figures [5] and [6l where the maximum lofting height reached, and 
the total time spent lofted, is shown for the full range of seismic simulations 
assuming that the slope with respect to the local surface gravity is 0° and the 
asteroids are non-rotating. Close to the micro-meteoroid impact the regolith 
is lofted on the p-wave and the lofting time and height decrease as a function 
of distance from the micro-meteoroid impact as the amplitude of the p-wave 
decreases. Then, at a certain distance, the lofting height and time suddenly 
increase; this corresponds to the distance at which the regolith starts to be 
lofted on the Rayleigh wave. This dip is particularly clear for the simulations 
with high attenuation (low Q). 

Note that the seismic simulations last 50 seconds and, for the lofting 
calculations, we take the assumption that at times greater than 50 seconds 
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the ground is stationary (this is a worst case estimate, especially for the low 
attenuation simulations). The ballistic regime is generally dominated by one 
long period of lofting (typically several 10s of seconds) that can be followed 
(and sometimes preceded) by several shorter jumps (see Fig. [7D- 

If the slope with respect to the local surface gravity is 0°, the horizon¬ 
tal motion of the regolith layer during the ballistic trajectory can also be 
estimated assuming that the regolith has a constant horizontal velocity de¬ 
termined by the horizontal (radial) velocity of the ground at the instant of 
ejection. However, the horizontal motion will not always be in the same di¬ 
rection (see example in Fig. ED- The total horizontal displacement for the 
regolith layer during the ballistic phase for all of the seismic simulations is 
given in figure d assuming that the slope with respect to the local surface 
gravity is 0° and the asteroids are non-rotating. For the low attenuation (Q = 
1000) seismic simulations, the maximum horizontal displacements due to bal¬ 
listic motion vary from ~6 cm at distances <10 m from the micro-meteoroid 
impact to ~0.01 mm at distances >300 m from the micro-meteoroid impact. 

If, however, the slope with respect to the local surface gravity is non¬ 
zero, the motion of the regolith layer during the ballistic trajectory will be 
strongly influenced by the direction of the local surface acceleration. As the 
component of the local surface gravity perpendicular to the slope will be 
reduced, the lofting times will be longer. Additionally, the component of the 
local surface gravity parallel to the slope will favour downslope horizontal 
motion (see Fig. ED- 

Our lofting calculations have shown that, during even a small micro¬ 
meteoroid impact (lOg at 6 km/s) the non-cohesive regolith layer on small 
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Figure 4: Examples of lofting behaviour, (a) The regolith layer is lofted by the p-wave 
and does not return into contact with the grjiipid until after the Rayleigh wave has passed, 
(b) The regolith layer is lofted by the p-wave and returns into contact with the ground as 
the Rayleigh wave is passing and is subsequently relofted. 


























Figure 5: Separation between ground and regolith during lofting - The maximum separa¬ 
tion achieved between the regolith layer and the ground as a result of lofting, as a function 
of distance (in m) to an impact of a lOg object at 6 km/s for homogeneous asteroid mod¬ 
els of following diameters: 100 m (magenta lines), 200 m (blue lines), 400 m (green lines) 
and 1 km (black lines); and following quality factors: 20 (top), 100 (middle) and 1000 
(bottom). The slope with respect to the local surface gravity is 0° and the asteroids are 
assumed to be non-rotating. 
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Distance from impact (m) 


Figure 6: Total lofting time - The total time the regolith spends lofted, as a function of 
distance (in m) to an impact of a lOg object at 6 km/s for homogeneous asteroid models of 
following diameters: 100 m (magenta lines), 200 m (blue lines), 400 m (green lines) and 1 
km (black lines); and following quality factors: 20 (top), 100 (middle) and 1000 (bottom). 
The slope with respect to the local surface gravity is 0° and the asteroids are assumed to 
be non-rotating. 
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Figure 7: Horizontal regolith motion during the ballistic phases - Top: the ground and 
regolith vertical displacements as a function of time showing the the ballistic periods. 
Middle: the radial velocity of the ground as a function of time. Bottom: The resulting 
horizontal displacement of the regolith during the ballistic periods. Note that the horizon¬ 
tal motion is not always in the same direction as it depends on the sign of the ground radial 
velocity at the moment of lofting. In this example the asteroid is 1 km in diameter, Q 
= 1000, the epicentral distance from the micro-meteoroid impact is 40° (700 m distance), 
the slope with respect to the local surface gravity is 0° and the asteroid is non-rotating. 
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Figure 8: Regolith radial displacement - The total radial displacement of the regolith 
during lofting , as a function of distance (in m) to an impact of a lOg object at 6 km/s for 
homogeneous asteroid models of following diameters: 100m (magenta lines), 200 m (blue 
lines), 400m (green lines) and 1 km (black lines); and following quality factors: 20 (top), 
100 (middle) and 1000 (bottom). It is assumed that the slope with respect to the local 
surface gravity is 0° and the asteroid is non-rotating. 
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R = 200 m, Q = 1000, distance from source = 70 m 



Figure 9: Influence of slope on regolith radial displacement - The total radial displacement 
of the regolith during lofting as a function of the slope with respect to the local surface 
gravity ( g ). The example given is for regolith at a distance of 70 m from a 
micro-meteoroid impact on the surface of a non-rotating asteroid of diameter 
400 m. The quality factor is 1000. The schematics explain the motion of the 
regolith during lofting. The initial velocity ( Vi ) is identical in all cases (with 
respect to the slope), and the subsequent trajectory is strongly influenced by 
the direction of the local surface acceleration, favouring downslope motion. 
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(sub-kilometric) asteroids is uplifted but not lost. The initial motion of 
the regolith layer will be dictated by the ground velocity at the moment of 
ejection. However, the following ballistic motion is strongly influenced by the 
direction of the ambient surface acceleration relative to the surface, favouring 
downslope regolith displacements. 

Similarly, there are other types of behaviour that may contribute to the 
radial motion. When the regolith layer is in contact with the ground and is 
subjected to horizontal motion coming from the ground (i.e., applied at the 
base of the regolith layer) downslope motion wi ll occur in a stop -start fashion, 


referred to as the stick-slip motion by 


Richardson et al 


(120051 ). H 


unde r stand the regolith motion the simple s l iding block analysis 


owever, to 


1965 


Richardson et al. 


2005 


Jibson et al. 


Newmarkj, 


2000|) is not sufficient because 


the seismic wavelengths in our simulations are likely to be of similar size, 
or perhaps even smaller than the depth of the regolith layer (see Fig. [Q 
Vp = 900m/s) . Aco ustic fluidisation of granular materials by the seismic 


waves (iMe l os h. 


19791 ) may also increase regolith motion in the same, downs¬ 


lope, direction. More detailed numerical simulations are required to study 
these types of regolith motion and any regolith flow that may occur as the 
regolith layer comes into contact with the ground at the end of a ballistic 
period. However, it is clear that over time, the cumulative effect of these 
mechanisms will result in accumulations of regolith material in surface ac¬ 
celeration potential lows or, “potential valleys”, as has been observed on 
several asteroids. On the surface of Itokawa, for example, the Muses Sea and 
Sagamihara regions coincide with the potential valleys and are generally ho¬ 
mogeneous, and relatively flat (slopes < 8°), consistent with a regolith layer 
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that has been allowed to seek out its minimum energy configuration after 


the formati on of the asteroid ({Miyam oto et ah 


Yano et al 


2006]) 


200 ?]; 


Fuiiwara et al. 


2006 


4. Implications of regolith concentration 


The non-cohesive regolith layer has a thermal inertia 3 to 4 orders of mag¬ 
nitude smaller than bare rocks forming the cohesive part (IDelbo’ et all 2007!). 
Our modelling predicts that, for a given source (here a micro-meteoroid im¬ 
pact), the induced surface accelerations tend to increase as the asteroid size 
decreases due to a decrease of the average distance to the source. Smaller 
asteroids will also have longer lofting times resulting in an increased mobil¬ 
ity of the loose regolith towards regions of low surface acceleration potential. 
However, as already pointed out above, nonlinearities in both the asteroid 
response to impacts and the uplift phenomena prevent a simple scaling with 
asteroid size. 

The high mobility of uncohesive regolith may expose an increasing num¬ 
ber of bare rocks of high thermal inertia thus helping to explain the ob¬ 
serve d trend o f in creasi ng average thermal inertia with decreasing asteroid 
size (IDelbo’ et ah . 120071); because the thermal inertia ratio between fine dust 
and bare rocks is so large, the surface average is strongly influenced by high 
thermal inertia regions. 

As a companion effect, the concentration of surface regolith in specific areas 
would create large lateral variations of thermal inertia from pol e to e q uator 


that may influence the Yarkovsky effect on these objects (IRubincam . 


Farinella et ah. 

1998; 

Bottke et ah, 

CD 

O 

O 

cm! 


1995 


2006 ). In the particular case of small fast 
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rotators, the regolith will concentrate in equatorial areas (surface accelera¬ 
tion potential lows), and consequently give this region a low surface thermal 
inertia. The sensitivity of such asteroids to the diurnal Yarkovsky effect may, 
therefore, be significant because the thermal energy will be re-emitted by the 
surface in a short delay comparable to the rotation period (Fig. UUT a)). More¬ 
over, the polar areas will be depleted in regolith and, consequently, will have 
a larger surface thermal inertia. Therefore, these objects may also be sensi¬ 
tive to the seasonal Yarkovsky effect, because the time necessary to re-emit 
the thermal energy is long, and comparable to the orbit period (Fig. [TOf bfh 
Consequently, even if the Yarkovsky effect depends on the pole orientation 
of the asteroid, smaller and faster rotators may have an increased orbital 
mobility due to surface thermal inertia variations induced by fine grained 
regolith concentration in equatorial areas. However, the significance of 
this effect as a function of asteroid size and rotation speed should 
be studied in more detail. 

In addition, an increase in regolith migration following micro-meteoroid im¬ 


pacts may enhance t 


l e effi ciency of regolith production via both impacts 


and thermal fatigue fiDelbo et all 2014) by regularly relocating the freshly 


produced regolith to the potential lows and, in turn, exposing more material 
to the surface thermal environment and to impacting meteoroids. 

The results presented here also have implications for lander missions on the 
surface of these small bodies, because a rough estimate of the regolith thick¬ 
ness is important for landing or penetration dynamics. In addition, as the 
lofting time computations presented in figure [6] do not depend on the lofted 
mass, these computations can be applied to any lander poorly coupled to 
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Thermal drag Thermal emission Spin direction Thermal Inertia 




Figure 10: Sketches of diurnal (a) and seasonal (b) Yarkovsky effects with latitudinal 
dependence of thermal ine rtia , (a) ob liq uity angle is 180° (retrograde), (b) obliquity 


angle is 90°. Adapted from 


Bottkeetal 


( 200611 . 


the asteroid surface. Therefore, the ground coupling capabilities and active 
source types of future seismic experiments are also affected by our results, 
which suggest that active sources may uplift uncoupled landers during a du¬ 
ration long enough to miss most of the seismic signal. 


5. Conclusion 

Assuming that asteroid interiors are cohesive enough to support the prop¬ 
agation of the seismic wave generated by micro-meteoroid impacts without 
irreversible deformations, continuum mechanics simulations of seismic wave 
propagation are presented for asteroids in the kilometric size range. Despite 
conservative hypotheses on the source amplitude and wave attenuation, we 
find that extensive seismic activity can be activated by the impact of even 
small micro-meteoroids. Therefore, seismic activity generated by impacts is 
much more frequent than assumed during previous studies. 
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Additionally, the high frequency surface accelerations generated by the 
direct seismic waves exceed the maximum surface gravity at a global scale for 
micro-meteoroid impactors (lOg at 6 krn/s). Therefore, during such events, 
any loose, uncohesive regolith at the surface of the asteroid is uplifted but not 
lost. The lofting times and displacements are estimated assuming a simplified 
ballistic behaviour as a first approximation and they present a non-linear 
behaviour depending on the seismic wave at the origin of the lofting process. 
During lofting the loose regolith is preferentially transported to regions of 
low surface acceleration potential. 

The concentration of fine grained regolith in specific areas has important 
implications for the thermal behaviour of the object. For example, the or¬ 
bital mobility of small fast rotators may be increased due to an increased 
sensitivity to both diurnal and seasonal Yarkovsky effects. We also present 
other implications for the deployment of landers or seismic experiments on 
the surface of these objects. 

Numerous possible implications of non-cohesive regolith uplift and concentra¬ 
tion induced by micro-meteoroid impacts on kilometric-sized asteroids have 
been discussed. However, this work can be continued in many directions: 
modelling of seismic wave scattering, numerical and experimental charac¬ 
terisation of loose regolith mechanical properties, including cohesive forces 


i n our model that cou l d play an importan 


(jHartzell and Scheeresl , 


2011 


Hartzell et ah 


role in the regolith behaviour 


20131) . modelling lateral varia¬ 


tions of thermal inertia and the consequences on the Yarkovsky effect, de¬ 
veloping the models to estimate the duration of the shaking following an 
impact and performing a detailed study of the timescales involved in global 
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redistribution of regolith (and crater erasure) given these predicted seismic 
surface accelerations and velocities. 
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